Melioidosis in Patients with COVID-19 Exposed to Contaminated Tap Water, Thailand, 2021.

In September 2021, a total of 25 patients diagnosed with COVID-19 developed acute melioidosis after (median 7 days) admission to a COVID-19 field hospital in Thailand. Eight nonpotable tap water samples and 6 soil samples were culture-positive for Burkholderia pseudomallei. Genomic analysis suggested contaminated tap water as the likely cause of illness.

Benchmarking CRISPR-BP34 for point-of-care melioidosis detection in low-income and middle-income countries: a molecular diagnostics study.

BACKGROUND: Melioidosis is a neglected but often fatal tropical disease. The disease has broad clinical manifestations, which makes diagnosis challenging and time consuming. To improve diagnosis, we aimed to evaluate the performance of the CRISPR-Cas12a system (CRISPR-BP34) to detect Burkholderia pseudomallei DNA across clinical specimens from patients suspected to have melioidosis. METHODS: We conducted a prospective, observational cohort study of adult patients (aged ≥18 years) with melioidosis at Sunpasitthiprasong Hospital, a tertiary care hospital in Thailand. Participants were eligible for inclusion if they had culture-confirmed B pseudomallei infection from any clinical samples. Data were collected from patient clinical records and follow-up telephone calls. Routine clinical samples (blood, urine, respiratory secretion, pus, and other body fluids) were collected for culture. We documented time taken for diagnosis, and mortality at day 28 of follow-up. We also performed CRISPR-BP34 detection on clinical specimens collected from 330 patients with suspected melioidosis and compared its performance with the current gold-standard culture-based method. Discordant results were validated by three independent qualitative PCR tests. This study is registered with the Thai Clinical Trial Registry, TCTR20190322003. FINDINGS: Between Oct 1, 2019, and Dec 31, 2022, 876 patients with culture-confirmed melioidosis were admitted or referred to Sunpasitthiprasong Hospital, 433 of whom were alive at diagnosis and were enrolled in this study. Median time from sample collection to diagnosis by culture was 4·0 days (IQR 3·0-5·0) among all patients with known survival status at day 28, which resulted in delayed treatment. 199 (23%) of 876 patients died before diagnosis and 114 (26%) of 433 patients in follow-up were treated, but died within 28 days of admission. To test the CRISPR-BP34 assay, we enrolled and collected clinical samples from 114 patients with melioidosis and 216 patients without melioidosis between May 26 and Dec 31, 2022. Application of CRISPR-BP34 reduced the median sample-to-diagnosis time to 1·1 days (IQR 0·7-1·5) for blood samples, 2·3 h (IQR 2·3-2·4) for urine, and 3·3 h (3·1-3·4) for respiratory secretion, pus, and other body fluids. The overall sensitivity of CRISPR-BP34 was 93·0% (106 of 114 samples [95% CI 86·6-96·9]) compared with 66·7% (76 of 114 samples [57·2-75·2]) for culture. The overall specificity of CRISPR-BP34 was 96·8% (209 of 216 samples [95% CI 93·4-98·7]), compared with 100% (216 of 216 samples [98·3-100·0]) for culture. INTERPRETATION: The sensitivity, specificity, speed, and window of clinical intervention offered by CRISPR-BP34 support its prospective use as a point-of-care diagnostic tool for melioidosis. Future development should be focused on scalability and cost reduction. FUNDING: Chiang Mai University Thailand and Wellcome Trust UK.

Genetic diversity, determinants, and dissemination of Burkholderia pseudomallei lineages implicated in melioidosis in northeast Thailand.

Melioidosis is an often-fatal neglected tropical disease caused by an environmental bacterium Burkholderia pseudomallei. However, our understanding of the disease-causing bacterial lineages, their dissemination, and adaptive mechanisms remains limited. To address this, we conducted a comprehensive genomic analysis of 1,391 B. pseudomallei isolates collected from nine hospitals in northeast Thailand between 2015 and 2018, and contemporaneous isolates from neighbouring countries, representing the most densely sampled collection to date. Our study identified three dominant lineages with unique gene sets enhancing bacterial fitness, indicating lineage-specific adaptation strategies. Crucially, recombination was found to drive lineage-specific gene flow. Transcriptome analyses of representative clinical isolates from each dominant lineage revealed heightened expression of lineage-specific genes in environmental versus infection conditions, notably under nutrient depletion, highlighting environmental persistence as a key factor in the success of dominant lineages. The study also revealed the role of environmental factors - slope of terrain, altitude, direction of rivers, and the northeast monsoons - in shaping B. pseudomallei geographical dispersal. Collectively, our findings highlight persistence in the environment as a pivotal element facilitating B. pseudomallei spread, and as a prelude to exposure and infection, thereby providing useful insights for informing melioidosis prevention and control strategies.

BurkHostGEN: a study protocol for evaluating variations in the Burkholderia pseudomallei and host genomes associated with melioidosis infection.

Background: Melioidosis is a frequently fatal disease caused by an environmental bacterium Burkholderia pseudomallei. The disease is prevalent in northeast Thailand, particularly among rice field farmers who are at risk of bacterial exposure through contact with contaminated soil and water. However, not all exposure results in disease, and infection can manifest diverse outcomes. We postulate that genetic factors, whether from the bacterium, the host or the combination of both, may influence disease outcomes. To address this hypothesis, we aim to collect, sequence, and analyse genetic data from melioidosis patients and controls, along with isolates of B. pseudomallei obtained from patients. Additionally, we will study the metagenomics of the household water supply for both patients and controls, including the presence of B. pseudomallei. Methods: BurkHostGEN is an ongoing observational study being conducted at Sunpasitthiprasong Hospital, Ubon Ratchathani, Thailand. We are obtaining consent from 600 melioidosis patients and 700 controls, spanning both sexes, to collect 1 mL of blood for host DNA analysis, 3 mL of blood for RNA analysis, as well as 5 L of household water supply for metagenomic analysis. Additionally, we are isolating B. pseudomallei from the melioidosis patients to obtain bacterial DNA. This comprehensive approach will allow us to identify B. pseudomallei and their paired host genetic factors associated with disease acquisition and severity. Ethical approvals have been obtained for BurkHostGEN. Host and bacterial genetic data will be uploaded to European Genome-Phenome Archive (EGA) and European Nucleotide Archive (ENA), respectively. Conclusions: BurkHostGEN holds the potential to discover bacterial and host genetic factors associated with melioidosis infection and severity of illness. It can also support various study designs, including biomarker validation, disease pathogenesis, and epidemiological analysis not only for melioidosis but also for other infectious diseases.

Phenotypic and genetic alterations of Burkholderia pseudomallei in patients during relapse and persistent infections.

The bacterium Burkholderia pseudomallei is the causative agent of melioidosis, a severe tropical disease associated with high mortality and relapse and persistent infections. Treatment of melioidosis requires prolonged antibiotic therapy; however, little is known about relapse and persistent infections, particularly the phenotypic and genetic alterations of B. pseudomallei in patients. In this study, we performed pulsed-field gel electrophoresis (PFGE) to compare the bacterial genotype between the initial isolate and the subsequent isolate from each of 23 suspected recurrent and persistent melioidosis patients in Northeast Thailand. We used whole-genome sequencing (WGS) to investigate multilocus sequence types and genetic alterations of within-host strain pairs. We also investigated the bacterial phenotypes associated with relapse and persistent infections, including multinucleated giant cell (MNGC) formation efficiency and intracellular multiplication. We first identified 13 (1.2%) relapse, 7 (0.7%) persistent, and 3 (0.3%) reinfection patients from 1,046 survivors. Each of the 20 within-host strain pairs from patients with relapse and persistent infections shared the same genotype, suggesting that the subsequent isolates arise from the infecting isolate. Logistic regression analysis of clinical data revealed regimen and duration of oral antibiotic therapies as risk factors associated with relapse and persistent infections. WGS analysis demonstrated 17 within-host genetic alteration events in 6 of 20 paired isolates, including a relatively large deletion and 16 single-nucleotide polymorphism (stocktickerSNP) mutations distributed across 12 genes. In 1 of 20 paired isolates, we observed significantly increased cell-to-cell fusion and intracellular replication in the second isolate compared with the initial isolate from a patient with persistent infection. WGS analysis suggested that a non-synonymous mutation in the tssB-5 gene, which encoded an essential component of the type VI secretion system, may be associated with the increased intracellular replication and MNGC formation efficiency of the second isolate of the patient. This information provides insights into genetic and phenotypic alterations in B. pseudomallei in human melioidosis, which may represent a bacterial strategy for persistent and relapse infections.

A Multiplexed Cas13-Based Assay with Point-of-Care Attributes for Simultaneous COVID-19 Diagnosis and Variant Surveillance.

Point-of-care (POC) nucleic acid detection technologies are poised to aid gold-standard technologies in controlling the COVID-19 pandemic, yet shortcomings in the capability to perform critically needed complex detection-such as multiplexed detection for viral variant surveillance-may limit their widespread adoption. Herein, we developed a robust multiplexed clustered regularly interspaced short palindromic repeats (CRISPR)-based detection using LwaCas13a and PsmCas13b to simultaneously diagnose severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and pinpoint the causative SARS-CoV-2 variant of concern (VOC)-including globally dominant VOCs Delta (B.1.617.2) and Omicron (B.1.1.529)-all the while maintaining high levels of accuracy upon the detection of multiple SARS-CoV-2 gene targets. The platform has several attributes suitable for POC use: premixed, freeze-dried reagents for easy use and storage; convenient direct-to-eye or smartphone-based readouts; and a one-pot variant of the multiplexed detection. To reduce reliance on proprietary reagents and enable sustainable use of such a technology in low- and middle-income countries, we locally produced and formulated our own recombinase polymerase amplification reaction and demonstrated its equivalent efficiency to commercial counterparts. Our tool-CRISPR-based detection for simultaneous COVID-19 diagnosis and variant surveillance that can be locally manufactured-may enable sustainable use of CRISPR diagnostics technologies for COVID-19 and other diseases in POC settings.

Antibiotic susceptibility of clinical Burkholderia pseudomallei isolates in northeast Thailand during 2015-2018 and the genomic characterization of ß-lactam-resistant isolates.

Melioidosis is an often fatal infection in tropical regions caused by an environmental bacterium, Burkholderia pseudomallei Current recommended melioidosis treatment requires intravenous ß-lactam antibiotics such as ceftazidime (CAZ), meropenem (MEM) or amoxicillin-clavulanic acid (AMC) and oral trimethoprim-sulfamethoxazole. Emerging antibiotic resistance could lead to therapy failure and high mortality. We performed a prospective multicentre study in northeast Thailand during 2015-2018 to evaluate antibiotic susceptibility and characterize ß-lactam resistance in clinical B. pseudomallei isolates. Collection of 1,317 B. pseudomallei isolates from patients with primary and relapse infections were evaluated for susceptibility to CAZ, imipenem (IPM), MEM and AMC. ß-lactam resistant isolates were confirmed by broth microdilution method and characterized by whole genome sequence analysis, penA expression and ß-lactamase activity. The resistant phenotype was verified via penA mutagenesis. All primary isolates were IPM-susceptible but we observed two CAZ-resistant and one CAZ-intermediate resistant isolates, two MEM-less susceptible isolates, one AMC-resistant and two AMC-intermediate resistant isolates. One of 13 relapse isolates was resistant to both CAZ and AMC. Two isolates were MEM-less susceptible. Strains DR10212A (primary) and DR50054E (relapse) were multi-drug resistant. Genomic and mutagenesis analyses supplemented with gene expression and ß-lactamase analyses demonstrated that CAZ-resistant phenotype was caused by PenA variants: P167S (N=2) and penA amplification (N=1). Despite the high mortality rate in melioidosis, our study revealed that B. pseudomallei isolates had a low frequency of ß-lactam resistance caused by penA alterations. Clinical data suggest that resistant variants may emerge in patients during antibiotic therapy and be associated with poor response to treatment.

Pneumococcal within-host diversity during colonization, transmission and treatment.

Characterizing the genetic diversity of pathogens within the host promises to greatly improve surveillance and reconstruction of transmission chains. For bacteria, it also informs our understanding of inter-strain competition and how this shapes the distribution of resistant and sensitive bacteria. Here we study the genetic diversity of Streptococcus pneumoniae within 468 infants and 145 of their mothers by deep sequencing whole pneumococcal populations from 3,761 longitudinal nasopharyngeal samples. We demonstrate that deep sequencing has unsurpassed sensitivity for detecting multiple colonization, doubling the rate at which highly invasive serotype 1 bacteria were detected in carriage compared with gold-standard methods. The greater resolution identified an elevated rate of transmission from mothers to their children in the first year of the child's life. Comprehensive treatment data demonstrated that infants were at an elevated risk of both the acquisition and persistent colonization of a multidrug-resistant bacterium following antimicrobial treatment. Some alleles were enriched after antimicrobial treatment, suggesting that they aided persistence, but generally purifying selection dominated within-host evolution. Rates of co-colonization imply that in the absence of treatment, susceptible lineages outcompeted resistant lineages within the host. These results demonstrate the many benefits of deep sequencing for the genomic surveillance of bacterial pathogens.

Highly specific and sensitive detection of Burkholderia pseudomallei genomic DNA by CRISPR-Cas12a.

Detection of Burkholderia pseudomallei, a causative bacterium for melioidosis, remains a challenging undertaking due to long assay time, laboratory requirements, and the lack of specificity and sensitivity of many current assays. In this study, we are presenting a novel method that circumvents those issues by utilizing CRISPR-Cas12a coupled with isothermal amplification to identify B. pseudomallei DNA from clinical isolates. Through in silico search for conserved CRISPR-Cas12a target sites, we engineered the CRISPR-Cas12a to contain a highly specific spacer to B. pseudomallei, named crBP34. The crBP34-based detection assay can detect as few as 40 copies of B. pseudomallei genomic DNA while discriminating against other tested common pathogens. When coupled with a lateral flow dipstick, the assay readout can be simply performed without the loss of sensitivity and does not require expensive equipment. This crBP34-based detection assay provides high sensitivity, specificity and simple detection method for B. pseudomallei DNA. Direct use of this assay on clinical samples may require further optimization as these samples are complexed with high level of human DNA.

Co-evolutionary Signals Identify Burkholderia pseudomallei Survival Strategies in a Hostile Environment.

The soil bacterium Burkholderia pseudomallei is the causative agent of melioidosis and a significant cause of human morbidity and mortality in many tropical and subtropical countries. The species notoriously survives harsh environmental conditions but the genetic architecture for these adaptations remains unclear. Here we employed a powerful combination of genome-wide epistasis and co-selection studies (2,011 genomes), condition-wide transcriptome analyses (82 diverse conditions), and a gene knockout assay to uncover signals of "co-selection"-that is a combination of genetic markers that have been repeatedly selected together through B. pseudomallei evolution. These enabled us to identify 13,061 mutation pairs under co-selection in distinct genes and noncoding RNA. Genes under co-selection displayed marked expression correlation when B. pseudomallei was subjected to physical stress conditions, highlighting the conditions as one of the major evolutionary driving forces for this bacterium. We identified a putative adhesin (BPSL1661) as a hub of co-selection signals, experimentally confirmed a BPSL1661 role under nutrient deprivation, and explored the functional basis of co-selection gene network surrounding BPSL1661 in facilitating the bacterial survival under nutrient depletion. Our findings suggest that nutrient-limited conditions have been the common selection pressure acting on this species, and allelic variation of BPSL1661 may have promoted B. pseudomallei survival during harsh environmental conditions by facilitating bacterial adherence to different surfaces, cells, or living hosts.

Corrigendum: An Evolutionary Arms Race Between Burkholderia pseudomallei and Host Immune System: What Do We Know?

[This corrects the article DOI: 10.3389/fmicb.2020.612568.].

Horizontal gene transfer rate is not the primary determinant of observed antibiotic resistance frequencies in Streptococcus pneumoniae.

The extent to which evolution is constrained by the rate at which horizontal gene transfer (HGT) allows DNA to move between genetic lineages is an open question, which we address in the context of antibiotic resistance in Streptococcus pneumoniae. We analyze microbiological, genomic, and epidemiological data from the largest-to-date sequenced pneumococcal carriage study in 955 infants from a refugee camp on the Thailand-Myanmar border. Using a unified framework, we simultaneously test prior hypotheses on rates of HGT and a key evolutionary covariate (duration of carriage) as determinants of resistance frequencies. We conclude that in this setting, there is little evidence of HGT playing a major role in determining resistance frequencies. Instead, observed resistance frequencies are best explained as the outcome of selection acting on a pool of variants, irrespective of the rate at which resistance determinants move between genetic lineages.

An Evolutionary Arms Race Between Burkholderia pseudomallei and Host Immune System: What Do We Know?

A better understanding of co-evolution between pathogens and hosts holds promise for better prevention and control strategies. This review will explore the interactions between Burkholderia pseudomallei, an environmental and opportunistic pathogen, and the human host immune system. B. pseudomallei causes "Melioidosis," a rapidly fatal tropical infectious disease predicted to affect 165,000 cases annually worldwide, of which 89,000 are fatal. Genetic heterogeneities were reported in both B. pseudomallei and human host population, some of which may, at least in part, contribute to inter-individual differences in disease susceptibility. Here, we review (i) a multi-host-pathogen characteristic of the interaction; (ii) selection pressures acting on B. pseudomallei and human genomes with the former being driven by bacterial adaptation across ranges of ecological niches while the latter are driven by human encounter of broad ranges of pathogens; (iii) the mechanisms that generate genetic diversity in bacterial and host population particularly in sequences encoding proteins functioning in host-pathogen interaction; (iv) reported genetic and structural variations of proteins or molecules observed in B. pseudomallei-human host interactions and their implications in infection outcomes. Together, these predict bacterial and host evolutionary trajectory which continues to generate genetic diversity in bacterium and operates host immune selection at the molecular level.

Genetic variation associated with infection and the environment in the accidental pathogen Burkholderia pseudomallei.

The environmental bacterium Burkholderia pseudomallei causes melioidosis, an important endemic human disease in tropical and sub-tropical countries. This bacterium occupies broad ecological niches including soil, contaminated water, single-cell microbes, plants and infection in a range of animal species. Here, we performed genome-wide association studies for genetic determinants of environmental and human adaptation using a combined dataset of 1,010 whole genome sequences of B. pseudomallei from Northeast Thailand and Australia, representing two major disease hotspots. With these data, we identified 47 genes from 26 distinct loci associated with clinical or environmental isolates from Thailand and replicated 12 genes in an independent Australian cohort. We next outlined the selective pressures on the genetic loci (dN/dS) and the frequency at which they had been gained or lost throughout their evolutionary history, reflecting the bacterial adaptability to a wide range of ecological niches. Finally, we highlighted loci likely implicated in human disease.

Genome-wide epistasis and co-selection study using mutual information.

Covariance-based discovery of polymorphisms under co-selective pressure or epistasis has received considerable recent attention in population genomics. Both statistical modeling of the population level covariation of alleles across the chromosome and model-free testing of dependencies between pairs of polymorphisms have been shown to successfully uncover patterns of selection in bacterial populations. Here we introduce a model-free method, SpydrPick, whose computational efficiency enables analysis at the scale of pan-genomes of many bacteria. SpydrPick incorporates an efficient correction for population structure, which adjusts for the phylogenetic signal in the data without requiring an explicit phylogenetic tree. We also introduce a new type of visualization of the results similar to the Manhattan plots used in genome-wide association studies, which enables rapid exploration of the identified signals of co-evolution. Simulations demonstrate the usefulness of our method and give some insight to when this type of analysis is most likely to be successful. Application of the method to large population genomic datasets of two major human pathogens, Streptococcus pneumoniae and Neisseria meningitidis, revealed both previously identified and novel putative targets of co-selection related to virulence and antibiotic resistance, highlighting the potential of this approach to drive molecular discoveries, even in the absence of phenotypic data.

Melioidosis in Thailand: Present and Future.

A recent modelling study estimated that there are 2800 deaths due to melioidosis in Thailand yearly. The Thailand Melioidosis Network (formed in 2012) has been working closely with the Ministry of Public Health (MoPH) to investigate and reduce the burden of this disease. Based on updated data, the incidence of melioidosis is still high in Northeast Thailand. More than 2000 culture-confirmed cases of melioidosis are diagnosed in general hospitals with microbiology laboratories in this region each year. The mortality rate is around 35%. Melioidosis is endemic throughout Thailand, but it is still not uncommon that microbiological facilities misidentify Burkholderia pseudomallei as a contaminant or another organism. Disease awareness is low, and people in rural areas neither wear boots nor boil water before drinking to protect themselves from acquiring B. pseudomallei. Previously, about 10 melioidosis deaths were formally reported to the National Notifiable Disease Surveillance System (Report 506) each year, thus limiting priority setting by the MoPH. In 2015, the formally reported number of melioidosis deaths rose to 112, solely because Sunpasithiprasong Hospital, Ubon Ratchathani province, reported its own data (n = 107). Melioidosis is truly an important cause of death in Thailand, and currently reported cases (Report 506) and cases diagnosed at research centers reflect the tip of the iceberg. Laboratory training and communication between clinicians and laboratory personnel are required to improve diagnosis and treatment of melioidosis countrywide. Implementation of rapid diagnostic tests, such as a lateral flow antigen detection assay, with high accuracy even in melioidosis-endemic countries such as Thailand, is critically needed. Reporting of all culture-confirmed melioidosis cases from every hospital with a microbiology laboratory, together with final outcome data, is mandated under the Communicable Diseases Act B.E.2558. By enforcing this legislation, the MoPH could raise the priority of this disease, and should consider implementing a campaign to raise awareness and melioidosis prevention countrywide.

Pneumococcal Capsule Synthesis Locus cps as Evolutionary Hotspot with Potential to Generate Novel Serotypes by Recombination.

Diversity of the polysaccharide capsule in Streptococcus pneumoniae-main surface antigen and the target of the currently used pneumococcal vaccines-constitutes a major obstacle in eliminating pneumococcal disease. Such diversity is genetically encoded by almost 100 variants of the capsule biosynthesis locus, cps. However, the evolutionary dynamics of the capsule remains not fully understood. Here, using genetic data from 4,519 bacterial isolates, we found cps to be an evolutionary hotspot with elevated substitution and recombination rates. These rates were a consequence of relaxed purifying selection and positive, diversifying selection acting at this locus, supporting the hypothesis that the capsule has an increased potential to generate novel diversity compared with the rest of the genome. Diversifying selection was particularly evident in the region of wzd/wze genes, which are known to regulate capsule expression and hence the bacterium's ability to cause disease. Using a novel, capsule-centered approach, we analyzed the evolutionary history of 12 major serogroups. Such analysis revealed their complex diversification scenarios, which were principally driven by recombination with other serogroups and other streptococci. Patterns of recombinational exchanges between serogroups could not be explained by serotype frequency alone, thus pointing to nonrandom associations between co-colonizing serotypes. Finally, we discovered a previously unobserved mosaic serotype 39X, which was confirmed to carry a viable and structurally novel capsule. Adding to previous discoveries of other mosaic capsules in densely sampled collections, these results emphasize the strong adaptive potential of the bacterium by its ability to generate novel antigenic diversity by recombination.

Interacting networks of resistance, virulence and core machinery genes identified by genome-wide epistasis analysis.

Recent advances in the scale and diversity of population genomic datasets for bacteria now provide the potential for genome-wide patterns of co-evolution to be studied at the resolution of individual bases. Here we describe a new statistical method, genomeDCA, which uses recent advances in computational structural biology to identify the polymorphic loci under the strongest co-evolutionary pressures. We apply genomeDCA to two large population data sets representing the major human pathogens Streptococcus pneumoniae (pneumococcus) and Streptococcus pyogenes (group A Streptococcus). For pneumococcus we identified 5,199 putative epistatic interactions between 1,936 sites. Over three-quarters of the links were between sites within the pbp2x, pbp1a and pbp2b genes, the sequences of which are critical in determining non-susceptibility to beta-lactam antibiotics. A network-based analysis found these genes were also coupled to that encoding dihydrofolate reductase, changes to which underlie trimethoprim resistance. Distinct from these antibiotic resistance genes, a large network component of 384 protein coding sequences encompassed many genes critical in basic cellular functions, while another distinct component included genes associated with virulence. The group A Streptococcus (GAS) data set population represents a clonal population with relatively little genetic variation and a high level of linkage disequilibrium across the genome. Despite this, we were able to pinpoint two RNA pseudouridine synthases, which were each strongly linked to a separate set of loci across the chromosome, representing biologically plausible targets of co-selection. The population genomic analysis method applied here identifies statistically significantly co-evolving locus pairs, potentially arising from fitness selection interdependence reflecting underlying protein-protein interactions, or genes whose product activities contribute to the same phenotype. This discovery approach greatly enhances the future potential of epistasis analysis for systems biology, and can complement genome-wide association studies as a means of formulating hypotheses for targeted experimental work.

Global and regional dissemination and evolution of Burkholderia pseudomallei.

The environmental bacterium Burkholderia pseudomallei causes an estimated 165,000 cases of human melioidosis per year worldwide and is also classified as a biothreat agent. We used whole genome sequences of 469 B. pseudomallei isolates from 30 countries collected over 79 years to explore its geographic transmission. Our data point to Australia as an early reservoir, with transmission to Southeast Asia followed by onward transmission to South Asia and East Asia. Repeated reintroductions were observed within the Malay Peninsula and between countries bordered by the Mekong River. Our data support an African origin of the Central and South American isolates with introduction of B. pseudomallei into the Americas between 1650 and 1850, providing a temporal link with the slave trade. We also identified geographically distinct genes/variants in Australasian or Southeast Asian isolates alone, with virulence-associated genes being among those over-represented. This provides a potential explanation for clinical manifestations of melioidosis that are geographically restricted.